ON FUNCTIONS OF JACOBI AND WEIERSTRASS (I) 



YU. V. BREZHNEV 



Abstract. We present some new results in theory of classical ^-functions of Jacobi and u- 
functions of Weierstrass: ordinary differential equations (dynamical systems) and series expan- 
sions. The paper is basically organized as a stream of new formulas and constitutes, in re-casted 
form, some part of author's own handbook for elliptic/modular functions. 



1. Introduction 

Thcta- functions of Jacobi and the basis of Weierstrassian functions a,Q,p, p' occur in numerous 
theories and applications. Since their arising, this field was a subject of intensive studies and 
towards the end of XIX century the state of the art in this field could be characterized about as 
follows. If some formula is not of particular interest then it has already been obtained. The majority 
of results were obtained in works of Jacobi and Weierstrass themselves and their contemporaries and 
followers (Hermite, Kiepert, Hurwitz, Halphen, Frobenius, Fricke, and others). There is extensive 
bibliography, including handbooks, on the theory of elliptic/modular and related functions. Most 
thorough works are four volumes by Tannery & Molk [7], two volume set by Halphen with a 
posthumous edition of the third volume [3] and, of course, Werke by Weierstrass [S] and Jacobi 
[Hj. A large number of significant examples can be found in "A Course of Modern Analysis" by 
Whittaker. 

In this work we would like to present some results concerning the ^-functions of Jacobi and a- 
functions of Weierstrass. Namely, power series expansions of the functions and differential equations 
to be obeyed by them. In a separate work we shall expound another way of construction of the 
Jacobi-Weierstrass theory different from Jacobi's approach of inversion of a holomorphic elliptic 
integral and Weierstrass's construction of elliptic functions as doubly-periodic and meromorphic 
ones. 

Series expansion of elliptic/modular functions is a widely exploited apparatus which, due to 
analytic properties, makes it possible to get exact results. It is suffice to mention the function 
series for various ^-quotients, the famous number-theoretic g-series, and their consequences like 
"Moonshine Conjecture" and McKay-Thompson series. 

Proposed dynamical systems and their solutions have an independent interest because solutions 
of physically significant differential equations, if these are expressed through Jacobi's 9- or modular 
functions, are consequences of the equations presented below. First of all we should mention the 
theory of integrable nonlinear equations and dynamical systems. Differential properties of the 
"modular part" of Jacobi's functions are more transcendental and diverse. They gave rise to nice 
applications over the last decade known as monopoles and dynamical systems of the Halphen- 
Hitchin type [JJ, Chazy-Picard-Fuchs equations considered in works by Harnad, McKay, Ohyama 
and others. Our interest in this part of the problem was initiated by a problem of analytic 
description of moduli space of algebraic curves and differential-geometric structures on it. This 
problem is completely open, so that the offered "differential technique" can be useful in those 
situations where some classes of algebraic curves cover elliptic tori. We shall not touch this topic 
and restrict ourselves only to "book-keeping" part of the apparatus. It is of interest in its own 
right. Applications suggest themselves from the formulas. 
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In sections 3-10 we shall be keeping the following rule. Unless one mentioned explicitly or if 
nothing has been said about formula then the result is new 1 . In some cases, even though a result 
is obvious, consequences can be nontrivial. In particular, this concerns the differential equations 
on ^-functions with respect to their first argument. Differential identities sometimes occur in old 
literature (see for example jjj), but the differential closure is needed. In most cases the derivation 
of formulas is self-suggested once a given formula(s) has been displayed. For this reason we omit 
proofs and restrict ourselves to comments to applications or elucidations some technical details 
which arc important for calculations. 

Content of section 2 is standard and presented here to fix notation. List of bibliography, modern 
and classical, reduced to a minimum since, even in a shortened form, it would require mentioning 
tens important works. 

2. Definitions and notations 

2.1. Functions of Jacobi. Four functions 9 1 2 3 4 and their equivalents under notation with char- 
acteristics are defined by the following series: 



'ii 



'10 



y 00 ■ "3 



6 x {x\t) = -i {-l) k ^ k+ ^ iT e^ k+1 ^ ix 

fe = -oo 
fe=-oo 

(x|r) = Y, 



gfc -KIT ^2k7TlX 



'01 



fe^fc 7rir g2k-xix 



We shall use the designation: 9k = 9k{x\r). Values of ^-functions under x = are called 
constants: dk = i?&(t) = #a;(0|t). For handy formatting formulas we shall use twofold notation for 
9- functions with characteristics: 

9[;}(x\r) = 9 a0 (x\r) = f] e ^+f) 2 -+ 2 <#+ §)(*+#) . 

We consider arbitrary integral characteristics and therefore the functions 9 a p are always equal to 
±#1.2,3.4- Shifts by ^-periods lead to shifts of characteristics: 

0[j\{x + z + f t\t) = 9[^:\{x\T)-e*' im{x+i ^ +IfT) , (n, m) = 0,±l,±2,... 

Twofold shifts by ^-periods yield the law of transformation of 0-function into itself: 

9 a p{x + n + m T \T) = 9 a0 (x\T)-e*^ m2T+2mx \-l) na ~ m P , (n, m) = 0,±1,±2, . . . . 

Value of any ^-function at any ^-period is a certain ^-constant with an exponential multiplier: 

0[?](i + f t\t) =^["/ + :](r)-e^ im( ^ + ^ r) , (n,m)=0,±l,±2,... 



^This means that the formula was not found in the literature, but does not mean that it is really new one. 
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2.2. Functions of Weierstrass. Weierstrassian notation a(x\co, ui'), p(x\u), u)'), ... and a(x; g 2 , #3), 
p(x;g 2 ,g 3 ), ... is the commonly accepted. Due to known relations of homogeneity for functions 
ct, C, p, p', the two parameters (w, oj') or (52753) can be replaced by one r = ^- which is called 
modulus of elliptic curves. We shall designate that functions as follows 

a(x\r) = a(x\l,r), ((x\t) = £(x\l, r) p(x|r) = p(x|l, r) p'(a;|r) = p'(x\l, r) . 

Weierstrassian invariants g 2 , g 3 and A = g 2 — 27gf are functions of periods or modulus and 
defined by known formulas of Weierstrass-Eisenstein. These series are entirely unsuited for numeric 
computations. Hurwitz, in his dissertation (1881), found a nice transition to function scries. Such 
series are most effective for computations: 



2 (r) = 207T 4 — + V — - „ , g 3 {r) 





77-function of Weierstrass is defined by the formula rj{r) = £(1|1,t) and the following formula is 
used for computations: 



where the numbers (a, 6, c, d) are integer and ad — be = 1. Modular transformations are used not 
only in theory but also in practice since the value of modulus strongly affects the convergence of 
the function series. Moving r into a fundamental domain of modular group (the process is easily 
automatized), one obtains the values of r with the minimally allowable imaginary part 3(t) = 
At such "the most worst" point the series converge very fast. 

Three functions of Weierstrass a x are defined by the following expressions: 

a x (x\u,uj')= e ^0f§, M^I/j) , where A = 1,2, 3. 

The cr-function of Weierstrass, as a function of (a;, g 2 , g 3 ), satisfies the linear differential equations 
obtained by Weierstrass: 



0(7 , da da 
ox dg 2 dg 3 

d 2 a da 2 2 da 1 2 

— — ^ — 12 — q 9 — Oo x a = 

dx 2 y3 dg 2 3 y2 dg 3 12 y2 



(1) 



These equations make it possible to write down the recursive relation for coefficients Ck(g 2 ,g 3 ) of 
the power series of the function a: 

<x;g 2 ,g s ) (■■■'■ ' C I.; ' ••• •'• ^ • ••• • (C = 1,C 1 =0). (2) 

Two such recurrences are known. One is due to Halphen: 

C k = - / Z>C k - l - 1 -{k-l){2k-l)g 2 C k - 2 , where % = -12 ff3 — - \ g\ — . (3) 

The second was obtained by Weierstrass: 

OO , 

1 \ — sr- a (9 2 \ m (o \ n * 4m+6 ™ +1 J A.o = 1 

cr{x;g 2 ,g 3 )- ^ A m,n[j) (2g 3 j (4m + 6n + 1)! > ( , (m < 0, n < 0) ' 

m,n— 

A m ,„ = y (n + 1) v4 m _ 2 ,„+i + 3 (m + 1) A m+ i in _i - i (2m + 3n - l)(4m + 6n- 1) A m _i ; „. 
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Recently, in connection with theory of Kleinian er-functions, yet another recurrence was obtained 
0. Among all the recurrences, Weierstrassian one is the least expendable. There are only multi- 
plications of integers in it and its comparison efficiency rapidly grows under the growing of order 
of the expansions. Weierstrass proves separately that the numbers A mM are integers. 

2.3. Function of Dedekind. Due to coincidence the standard designations for Weierstrassian 
function r\(r) and Dedekind's one, we shall use, for the last one, the sign rj(r): 



rj{r) = e 12 



''IK* 



„2fc7I 



(3fe +fe)7ri 



(Euler (1748)) 



There is a differential relation between the two functions. It is expressed by the formula i — = — ri. 

<r) dr n 

3. Differentiation of functions a, C, p, p' 

Differentiations of Weierstrassian functions with respect to invariants are known |3J |7j. Rules 
of transformations between differentiations with respect to {g 2 ,g^) & n d (uj,lu') are also known 
(Frobenius-Stickelberger (1882)) ^j. Applying them, and assuming J > 0, one obtains 



da 

doj 


'-{^(p-e-^g^ 


da 
du>' 




(K = _ 

du> 


-i {u'(p' + 2(p-±g 2 


K = 

doj' 


- {^(p' + 2Cp- \ g 2 


dp 
duj 


^ {^'(2p 2 + Cp'-i. 


dp _ 
duj' 


~ [ui(2p 2 + Cp' -\t 


— \uj' 

7T I 


(6pp' + 12(p 2 ~g 2 C) - 



duo 

= _i{ w (6pp' + 12Cp 2 -5 2 C) --n{&p' + \2xp 2 -g 2 x)} 
Setting here (uj = 1, u>' = r) one obtains the dynamical system with a parameter x: 
da 



dr 

dr 

dp 
dr 



dr 



^{p-( 2 + 2 v (xC~l)~^g 2 x 2 }a 
- | p' + 27 1 ( + 2p{C-XT]) ~^52^} 
-2±- {2 p 2 + p'(C - xr,) - 2 77 p - i ,g 2 } 
-6^{p'(p-r7)+ (2p 2 - ift^fc-zjj)} 
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The function a(x\r) satisfies one differential equation 

d 2 a da . da / 1 2 \ 

— -2x V -~m- + (2 V +-g 2 x )a = 0, 

which is an analogue of heat equation for the function 9{x\t). All the differential equations in 
variables x,u),u}',t contain no the quantity g 3 . Depending on representation (g 2 , g 3 ), (uj,uj') or r, 
it is either algebraic integral of these equations g 3 = 4 p 3 — g 2 p — p' 2 or determines a fix algebraic 
relation between the variables. Each of the functions a, p, p'{x\r) satisfies ordinary differential 
equation of 4-th order with variable coefficients g 2 3 (t), t](t) in the variable r. These equations 
have defied simplifications. By virtue of linear connection between function a and 9, the last one 
also satisfies ordinary differential equations. 



4. Power series for ct-functions of Weierstrass 
Weierstrassian recurrence has the following graphical interpretation: 



n 





-2, 

^-v. (rn, n) 
m.. .^^i 




(m — 1 , n) 


\. (m + 1 , n — 1) 











o m 



The explicitly grouped cx-series © has the form ([n] denotes an integral part of the number) 



CO (k/2 



^{x;g 2 , 



r 2fe+l 



(2ft + 1)! 



».*>)=£ 1 E 2 2fc - 5 ^ 3 ^ fc „ 2 , . 9 r- fc 53 fe - 2u 

fc = U=[fc/3] 

Denote e A = p(wjcj,a;'). Then the functions er A satisfy the equations of Halphen: 



da. da. da. 

x^±-2e x ^~4g 2 -± = 

dg 2 



dx 



de x 



.&-(<* , -S«0&-^-**>&+(*+£*-'H 

Herefrom one obtains analogue of the recurrence (@J : 

OO ( k/2 



..2/.' 



(2ft)! ' 



(4) 



k=0 < "=° 

Q3 ro ,„ = 24 (n + 1) <B m _ 3)n+1 + (4 to - 12 »- 5) 58 m _i,„- 

-|(m + l)<B ro+ i,„_i - i(TO + 2n- l)(2m + 4n-3) Bm,n-i- 

All the coefficients 58 m)Tl are integrals, and 5B 0j0 = 1 and 58 m , « = under (to, n) < 0. Weierstrass 
wrote out recurrences for the functions = e^ exX a x in representation (e A , e A = 3e A 2 — jg 2 )- 
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One can build a universal series for the functions a, o~ x . All the functions satisfy the differential 
equations 

(5) 

0-( 4 * 2 -!*)f; -i2(W-**)g + ( e * + ^)3 = o 

where case S = o~ x corresponds to e — 1, and H = a corresponds to e = and arbitrary e A . The 
quantity e is a parity of ^-characteristic. The universal series and integer recurrence acquire the 
following form: 

00 r fc / 2 ^ 9fc-n- P 



k=o {»=° 



(2jfe + l- e)! ' 



!8£> n = 24 (n + 1) <CU, n+1 + (4m - 12 n - 4 - e) QC-i.n" ' ° 1 



(m+l)%% +hn _ 1 -Um + 2n-l)(2m + 4n-l-2e)<8% n _ 1 



l 



Transition between pairs (52133)' (^jS^) an d ( e A,e t ,) is one-to-one therefore the universal recur- 
rence can be written for any of these representations. In the last case it is symmetrical. The 
representation 10 is a five-term one, in contrast to four-term recurrence of Weierstrass A mj „, but 
the representation (<? 2 i33) 1& n °t possible for the functions a x . 

5. Power series for ^-functions of Jacobi 

Jacobi made an attempt to obtain such series as early as before appearance of his Fundamenta Nova 
[HJ I: 259-260]. These series are the series with polynomial coefficients in t](t) and 0(r)-constants. 
This fact follows from the formulas 

9£x\t) =irfj 3 (T)-e- 2riiT)x2 a(2x\T) , x (x\t) = ^•e^^o^zlr) . 

^-constants make it possible to rewrite Weierstrassian representations g 2 , <73,e A into ^-constant 
ones. Formulas for branch-points e k through the ^-constants are well known. In turn, ^-constants 
are connected via Jacobi identity $3 = d^ + ^X- All this enables one to convert the representations 
choosing arbitrary pair. We shall speak (a, (3) -representation 1 if formulas are written through the 
constants (0 QO , O /?) under (a,/?) ^ (0,0). Operator © has Weierstrassian representation 



S = -12( 4eA 3 - ff2eA )^--(4e A 2 -|g 2 ) 



= -(2ej + 2e,e x - e 2 )^- + J(2e 2 + 2e A e tl - e t 2 )^- , 
3 v ' de x 3 v ' cte M 

and also z9-constant one. For example (0 2 , ^^-representation 

Let us denote (a) = (—1)". General (a, /^-representation of Halphen's operator has the form 



© = y ((a) V% + 2 </3> tff O ^1 - y ((/3> t?So + 2 <a> 1?*o O • 
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Let e is defined by a parity of characteristic of the function Q/3 (a;|r): 

(a/3) + 1 | £ = under 9 aff = ± 6 1 



e = 



2 1 e = 1 under 9 al3 = ± 9 2 3 4 J 

The equation © for the function 2 = (cr, cr A ) has the form 

where 2) is determined by formula JJJ). 

5.1. Function 6^x\t). Expansion of function 

0^1 r) = Y, C k {r)-x 2k+1 =2^{«-2^ + (2^-^ + ^ + ^). B » + ...} (8) 

fe=0 

has the following analytic representation: 



) 1 (a ! |T)= 2tt^ 



fc=0 



(47ri) fc d fc ?7 3 
(2fc + 1)! dr fc 



(9) 



(-i)=© s ,„_ s .< s <- s) 



fc=0 <. "=° 

where polynomials Af u ("&), depending on combinations of ^-constants, have the form 

I *-*v — a, a " ' 

= J] 

Integer recurrence ® TOl „ (<S ,o = 1) looks as follows: 

S„ = 4(n-2m-l) <S m ,n-i -4(m-2n- 1) <8 m _i, n - 

- 2(m + n- l)(2m + 2n - l)(0 m - 2 ,„ + B -i,„-i + 8 m ,„-j) 

The recurrence © m .„ is antisymmetric: <S TO ,n = (—1)'™+" (S nm . One may write out representation 
of the type Q\(x\t) = X] C mnp c/™ g% if x k through the recurrence of Weierstrass A m<n , but <8 m ,n 
is more effective than A m n since polynomials are already grouped together in ^-constants. 

Odd derivatives 6*^ 2fc+1) (0|r) = $f h+1 ' (t), i.e. expressions in front of x 2k+1 in ©, generate 
polynomials (JSJ) in (77, 1?) which are exactly integrable k times in r. 

5.2. Functions #2,3,4(;e|t). Expansions of functions 6[} = ±#2,3,4 of the form 

00 2 

= E ^ M ■ a2 * = - { 2 »? + t ( </3> ^ V]4 " <Q> ^ [ A]4 ) } * 2 + ' ' ' (10) 

fc=0 

have the following analytic representation: 

The equation (JjjJ on functions S = (cr, tr A ) in (a, /^-representation has the form 
d 2 E 



dx 2 



S E + + ^ (02, + (a + /3) tff < + < 3 )^} S 
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where e^ s ($) correspond (independently of representation (a, j3)) to chosen function a or <r A : 

e, s m = g o> * m 4 - (7> ^A] 4 ) , { I- ; *° = m . (12) 

Representation is called symmetrical if (7,(5) = (a, (3). 

In symmetrical representation the series 111 111 has the following grouped form: 

DO f k \ 

fl[?](x|r) = <?[?](r) (-2)' E L'ImIV ^^^ (13) 
fc=o I -=° J 

with the universal integer recurrence (& a f } = 1): 



(«.» ,0f o l4« 9 r Q - 1 i4( ! y-s) 



s=0 

= (a>(4n-8m-3)©^ : ; ) _ 1 -< / 3>(4m-8n-3)C-l,n- 

- 2 (m + n - 1) (2 m + 2 n - 3) (©^> _ „ + (a + /3> 6^ „_ 1 + 0£;»_ 2 ) . 

Symmetry of the recurrence 25^'^ with respect to permutation of indices is determined by the 
properties: 

^n,m \ *■/ ^m,n ' 

rft(/3,o) _ f_1\(i»+n)(»+S) (RCo.fl) 

w m,n V / v m, n ' 

This means that there are only two recurrence ©^fn under a = {1,0} and (3 = 0: 
= <a>(4n-8m-3)0^„_ 1 -(4m-8n-3)©^ , _ 1 ,„- 

- 2 (m + n - 1) (2 m + 2 n - 3) (0^_ 2> „ + (a) &ZU n-x + <8£»_a) 
with property <S£f m = (-l)("»+"><«+» ®£ n . 

U5<°) = ('_i)(™+»)*(<>) <ft(i) = (Kt 1 ) 

^77,771 V / ^m,n > 7i, m ^777,77 • 

Developments (j5J) and (|13|l differ from each other only in a multiplier and shape of the recurrence. 
They can be unified into one (by introducing the parity e) but the quantity (a) is left and (m, re- 
entries of matrices <3 {/3,a) differ only in a sign. 

Even derivatives 9^ (0\t) = ^^'(r), i.e. expressions in front of x 2h in l|13ll . generate polyno- 
mials (|1U|I in (77,1?) which are integrable k times in r. Their integrability is a consequence of one 
dynamical system considered in sect. 7. 

Making use of these expansions one can build other ones about points x = {±A, ±5}. Described 
series will be, up to obvious modifications, changed into one another. 

6. Dynamical systems on ^-functions 
6.1. Differential equations in x. The five functions 

09 

01,2,3,4 and 6 1 = 



ON FUNCTIONS OF JACOBI AND WEIERSTRASS (I) 

satisfy closed ordinary autonomous differential equations: 

^ - ^-6 -7T9? 2 .^i 96»i_6»i 2 2 2 2 61 r ^ . •» 

<9#3 _ ,y l n2 u :r : ""1 _ /y 

<9# 



4 _ u l /) _„q2 



3^4 




h 


dx 


2#4 


de 1 


*1 


dx 


2^3 




*1 





-dx- = r 1 ^- n 

These equations are equivalent to relations between Weierstrassian functions (a, (, p, p')(x\r). Gen- 
eral form of differentiation of the functions #1,2,3,4 is as follows = 0): 

d0 k 0' 2 v e il 8fc-28 10fc-28 n , ,x 

It immediately follows that the following identities are obeyed 2 

Ql 01 Q 0^ 

-— ' — -p- = 7rt? 2 • 1 sign(n - m) , where (fc 7^ n, n ^ m, fc 7^ m, fc, n, m = 2, 3, 4) , 

and also well-known polynomial identities: 

# + = + tf 2 # 2 =tf 2 # 2 + tf 2 # 4 2 , sign(n-m).,9 2 2 = C^-^€- (14) 

6.2. Differential equations in r. The functions #1,2,3,4, 0[ satisfy closed non-autonomous ordi- 
nary differential equations in r: 



96*3 -i #i 2 i 2 a a Q'l Tri „ 2 ,a2 ^ 2 



tf§-|ff 8 + {>+gW + ^}-«8 • (15) 

General form of r-differentiation of the functions #1,2,3,4 is as follows: 

+ + where ,= |^f, , = i2£f>, . = (1.2.3,4). 



These fundamental relations appear implicitly in works by Jacobi but have not got into comprehensive handbook 
for elliptic functions compiled by Schwarz on the basis of Weierstrass's lectures. 
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7. Differential equations on 0, ^-constants 
Differential closure is provided by the constants $2,3,4 and 77: 

} 3 f - ;■> + 5 W - *8 ■ S = M 2 " 2 - £w + •>» + •>«)}■ 

Well known differential equations on logarithms of ratios 2 : $3 : $4 and also known dynamical 
system of Halphen and its numerous varieties are widely encountered in the modern literature. 
These are generated by the system of equations 



52, 



which, in implicit form, was written out by Weierstrass and Halphen used it I: p. 331] in order 
to get his famous symmetrical version of the system. Ramanujan obtained its equivalent for some 
number-theoretic g-series. See for example § 1] and additional information in this work. It is 
lesser known, to all appearances, that another version of the system was written out by Jacobi 
in connection with power series of ^-functions. These results were published by Borchardt on the 
basis of papers kept after Jacobi II: 383-398]. 

Jacobi obtained a dynamical system on four functions (A, B, a, b) which are, in our notation, 
rational functions of the 0, //-constants 

A = 02, B — — — —4r ~\~ — 2 -g 4 , = 4 1-2-1 , 6 = 2- 2 4 



3 $1 ' v $17' $3 

and showed that the dynamical system can be useful for obtaining power series of 0-functions. He 
also presented two (canonical) transformations of the variables (A, B, a, b) retaining the shape of 
the equations. Jacobi noticed that the developments become simple and have a recursive form 
when extracted an exponential multiplier q~2 ABx . This is, in fact, Halphen's recurrence © for 
the cr-function of Weierstrass. Transition between variables in mentioned systems is not one-to-one 
but always algebraic. This is because all of them are consequences of the equations (jlfij) since 
the variables, appeared in these systems, are rationally expressible through the 0, 77-constants. A 
related dynamical system arose in a work of Jacobi earlier, when he derived his known equation 
on ^-constants HO II: p. 176] 

(d 2 d TTT - 15 d d T d TT + 30 0*) 2 + 32 (tf d TT - 3 2 ) 3 = -TT 2 d w ($ t? tt - 3 t? 2 ) 2 . 

Logarithmic derivatives of 0-constants, in turn, satisfy compact differential equation 



(/ T - 2 f) f TTT - f TT + 16 ff TT + 4 f T (f T - 6 f) = 0, / = |- In 2i3 , 4 (r) 



and Dedekind's function A = In 77(1") satisfies the differential equation of the 3-rd order 
{ A rrr - 12 A rT A r + 16 A 3 T } 2 + 32 { A TT - 2 A 2 } 3 = A n e fj 2A . 
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8. Modular transformations of functions 6\ and rj 



The modular transformation is necessary for computations of 9, 77- functions. General PSL^Z) 
transformation for the function 8 1 is closed in itself (in contrast to 9 2 3 4 ): 



ZI±*) =e g (a,c,d)-VcT + d e^+X O x {x\t), e e (a,c,d) = v 7 ! 



cr -f d 



P>7ri 



e g (a,c,d) = e 



f(2c-3)-i 



^ S ign(-d) + i ^ k[ik] 



under normalization c > 



O x {x\t + N) = c^-O^t) 

When discussing the modular properties of ^-functions one is usually pointed out "where Sq is 
an eighth root of unity" but complete algorithm of its computation leaves aside. Meantime, as it 
follows from the shape of ^-series, it is seen that the series, which can be transformed into a "hyper- 
convergent" form, can turn out to be incomputable. It is known that the multiplier £q((i, C, d) can 
be expressed via Jacobi's symbol (|) (Hermite). It has independent rules of computations. These 
formulas are displayed here for completeness of computations and contain some simplifications of 
Dedekind's sums (see about them): 



ar + b 
ct + d 



= C 



a=f-| (2c-3)-i-^isign(-d) + I J2 k[ik] 



y/cr + d rj{r) 



rj(r + N) = e SJV J7(r) 



9. ^-FUNCTIONS WITH CHARACTERISTICS 

The preceding results, making use of notation with characteristics, allow us to unify formulas and 
can be a subject for further generalizations to higher genera. 



9.1. (a, /^-representations. Any object, symmetrical in ^-constants, can be rewritten in (a, /?)- 
representation. For example representation of branch points through the ^-constants (112(1 . $ Q|3 - 
representations for quantities g 2 3 have the following form: 



.92 ( 



(t)= ^ {#lo + (a + i3)#i + , 



(a,/3)^(0,0) 



9 3 (r) = ^ {2 {/3> 0% - 3 di (</3) <, - <a> < ) - 2 (a) j . 

Identity $3 = i?| + d\ and Jacobi's formula = 27tt/ 3 , in (a, /^-representation, have the form 



^] 4 =fw^„f + (a>^ <a/3> + 1 



(17) 



^(r) = i' 3 (1 - <a/3))-7ri^ 3 (r) , where ^(t) = 0^(O\t) 



In general case, when the ^'-constant is a value of the derivative of ^-function at some ^-period, 
the formula of Jacobi (|17f) is generalized to the following expression: 



n , m 



T =(-i)« 



•7rie 



2 r |./3+„ 



l-( a + m>^ +n )-^-m^]} . (18) 
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Algebraic and differential closeness of Jacobi's functions with integral characteristics lead to the 
fact that all the derivatives 9^ i3 {x\t), with arguments shifted by arbitrary ^-periods, are expressible 
in terms of the function 9' 1 (x\t) and functions 9 afj (x\T): 

«•(• + 1 + 7 +) = H)— - ™) *[S*M - 

— (a + m) ■■Ktt[ p+ „\ 



Under x = need to use the formula (fT%|) . 



9.2. Modular transformations. The modular transformations served as the key for Hermite 
when he obtained his celebrated solution of the quintic equation x 5 — x — A = 3 . Hermite wrote 
out the multiplier £ g (a,c,d) in a form of Gaussian sum of exponents. We propose the exponent 
of a sum what is simpler. General modular transformation r — > f^rpj for the functions 9 has the 
form 

f ^if— ^) =e e (a,c,d)e^ 2a{bcp - d+1) - cP{ ^ 

\CT + d CT + dj 



3Ti| fif-|(2c-3)-|-£=i S ign(-d) + i £ fe[|fe]j 
£g(a, c, d) = G fe=1 , under normalization c > 

t 0M(z|r + JV) = e^^^^^l^lr) 
where, at given modular transformation, characteristics are computed by the formulas: 

oi = da — c{3 ( a — aa' + cj3' 

/3' = -ba + a(3 \f3=ba' + d/3' 

Under x — the transformation turns into the transformation of the ^-constants 

^ ?'-d n = £ e {a,c,d e • V c r + d i9 (r) 



3Ti|&|-#(2c-3)-i-£ii S ign(-d) + i £ fc[e fe ]j 
£g(a, c, d) = e , under normalization c > 

[ l ?[7 1 ](r + AT)=e^ iV(a2 - 1) .^](T) 

If characteristics (a, /?) = {0, 1} mod 2 then the formulas are closed: (a', (3 1 ) = {0, 1} mod 2. Ratio 
of any two §, 9- functions contains no the multiplier £g(a, c, d). Hermite used this fact to build his 
functions ip, ip, x( T ) an( i transformations between them. 



''Solution of this equation is expressed, in fact, in terms of i9-constants and was obtained by Hermite in 1858 
(see last formula on p. 10 in Curiously that it contains an erroneous sign which is repeated, to the best of our 
knowledge, everywhere the solution is reproduced. 
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9.3. Multiplication theorems. Let n be arbitrary complex number. Then the functions 9 satisfy 
the complex multiplication theorems determined by the following recurrences : 



i{nx) = °J n ^ *2 (*) - Ol(n lX ) 01 (*) 



6» 1 (2a;) =20 1 (x) 



(fi) 0{f(n lX ) 0{f(x) + (a) e[f(n lX ) 0{f(x) 



-n 2 



)Xl]((n-2)x) 



Ux)e 3 (x)e i (x) 



n 1 = n — 1 



This is not the only representation. Multiplication for 9[ is derived by taking a derivative. If n 
is integer then the formulas are closed. In particular, multiplications for the functions #2, 3, 4 are 
closed (doublings were written out by Jacobi): 

el(n lX ) 9l{x) -0l(n lX ) 9\{x) 



9 2 (nx) 
9 3 (nx) 
9 4 (nx) 



$l-6 2 {{n-2)x) 

el{n lX ) 9l(x)+9l(n 1 x) 9j(x) 
#l-e 3 ((n-2)x) 

ej(n 1 x)9l{x)-ej{n 1 x)9l(x) 
$l-6 A {{n-2)x) 



There exist the polynomial non-recursive multiplications but these are multiplications not over the 
constants z? Q/3 (r) . This follows from Weierstrassian formulas for the functions er, a x (nx) . 

9.4. Differential equations. 

9.4.1. 9 -functions. 9- functions with arbitrary integer characteristics (a,/3), as functions of the two 
variables x and r, satisfy splitted and closed systems of ordinary differential equations over the 
differential field of ^-constants Cd($,rj): 



89 



dx 



6i [ 



d9[ 
dx 



(a) TT1 



^ - am n "I 1 



,2 g[V]fl[A 



(19) 



{ 4 r ? + ^(^; 4 + ^; 4 )}-c 



+ {>+5W: 4 +€ 4 )}^[;i • (20) 

t: §f + 3 { t^ 2 ^ • % + > + S era 4 + ^ 4 ) } * - 
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These differentiations are not quite symmetrical. In all likelihood, the separability of variables in 
differential properties of the 0's is not accidental. Especially in those cases, where exact solutions 
of integrable nonlinear partial differential equations are expressible through the ^-functions. All 
such solutions, essentially non-stationary and multi-phase, are consequences of integrability of the 
only system of ordinary differential equations on the five Jacobi's functions 9, 9[. They satisfy the 
heat equation but it is a partial differential equation. 

General solution of dynamical system 1)19(1 is given by the formulas 

e af3 = a ■ 6 a0 (x + b\r) e cx , 6[ = a ■ {9[(x + b\r) - c9 n (x + b\r)} e cx 
with arbitrary functions a, 6, c(r). The system (|20|) has a solution of the form 

9 aP = a- e a0 {b\r) , 6[ = -a' ■ 9 n (b\r) + a ■ ^(6|r) , 
with arbitrary functions a, b(x). All the solutions should be supplemented by the relations i|14[l. 
9.4.2. d-constants. (a, /^-representation for the system of equations l|16[) has the following form 

dm 



dr 



''" kW'^i + ^ + {a + p) ^ } *- (a ' /3) ^ (0 ' 0) 



(21) 



dr 

Its general solution, with three arbitrary constants (a : b : c : d), is given by the formulas 

'ar + b\ i _ , JN _2 „( ar + b \ , 1 n i° 

2cr + d 

and the first identity of Jacobi l|17jl. Right hand sides in these formulas are the 77-series. Unless 
one assumes that the quantities 77,1? in equations l|2()JI being the t](t), z?(r)-constants, then the 
requirement of differential closeness of the field of coefficients of the system, due to heat equation, 
leads to the equations (|21|) . One may also view the equations l|21|) as the integrability condition 
(compatibility) of the two systems of equations l|19() and (|20|l . 

10. Conclusion 

Dynamical systems of Jacobi-Halphen, their generalizations, known as SU (2)-invariant self-dual 
Einstein's equations [S], modular solutions of equations of Painlevc P VI (Picard, Okamoto, Manin, 
and others) etc are also consequences of splittability of the equations on ^-functions. Thus, rules of 
differentiations (|19II21|) generate mentioned dynamical systems and their solutions and the quan- 
tities 9, 77, i9 are the uniformizing variables for them. Higher derivatives of the ^-functions or 
^-constants are again the 77, 1), 9, 6^-functions in form of rational polynomials of them. Using these 
facts one can find new solutions or simplify known ones. For example solution of equation P V i 
obtained by Hitchin in the work 0. 

Modular transformations of sect. 9.2 generate automorphisms of the systems H19H21|I . For the 
7?, ^-functions these are merely linear transformations. For other dynamical systems, like Jacobi's 
equations on (A, B,a,b), this yields non-obvious transformations of dynamical variables and sys- 
tems into themselves. For the modular solutions of equations of the Painleve type the transforma- 
tions become lesser obvious and cease to be canonical and change shape of the equations. 

In conclusion we note that it is of interest to generalize preceding results to 9-functions of 
higher genera and their constants. In some particular cases this problem might be tested in the 
following manner. Let us consider a situation when the Jacobian variety admits a decomposition 
of two-dimensional O-function to functions of Jacobi. For example 

e(S 1 1 l) = \ (U*i\t) + U*i\t)) 3 (z 2 \fx) + \ (9 3 ( Zi \t) - 9 4 ( Zi \t)) 4 (z 2 | m ) . 



ON FUNCTIONS OF JACOBI AND WEIERSTRASS (I) 



15 



Corresponding ten i9 [^-constants are expressible through the quantities 

l?2,3,4( T ); #2,3,4^), ^ 1,2,3,4(3 | t ) . 01,2,3,4(11*0 ■ 



Making use of the relations 



2 #3(3) = ^3 + , 2 = 4 0§ - ^ , 

1 $ 2 1?, $ 4 

#4(3) = ^(3) ' Hi) = 9 flTTTflSTTT 




(these are obtained from the multiplication theorems) one may analyze the system of derivatives 
of the t9-constants with respect to moduli (r, /i). In particular, to check its splittability in r 
and [i. This fact could be also checked, due to the heat equation, by investigating closeness of 
differentiations of sixteen functions Q\%}{z 1 ,Z2) with respect to arguments (zi, z 2 ). As it is seen 
from this example, the closure does exist and it being algebraic at the least. Whether it will be 
polynomial and, if no, how to do this? 
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